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融合 差分 变异 和 切线 飞行 的 天 鹰 优化 器 C 
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(上 海 工 程 技术 大 学 管理 学 院 ， 上 海 201620) 


摘 针对 天 应 优化 器 (aquila optimizer，AO) 虽 然 拥 有 强大 的 全 局 勘探 能 力 ， 但 局 部 开发 能 力 不 足 的 问题 ， 提 出 

融合 差分 变异 和 切线 飞行 的 天 鹰 优化 器 (differential evolution mutation and tangent flight aquila optimizer, DEtanAO). 

首先 ， 根 据 差分 进化 算法 中 的 变异 操作 能 使 算法 具有 较 强 的 开发 能 力 弥 补 了 AO 算法 的 不 足 ， 然 后 ， 利 用 切线 搜索 

算法 中 切线 飞行 策略 具有 较 强 的 探索 搜索 空间 的 能 力 并 能 使 算法 跳出 局 部 最 优 解 的 优势 ， 用 其 替换 了 AO 算法 中 的 

t 这 两 种 策略 的 结合 有 效 地 平衡 了 DEtanAO 算法 的 勘探 和 开发 阶段 。 最 后 ， 为 验证 DEtanAO 算法 的 优化 
生 能 ， 在 12 个 标准 基准 函数 、 高 维 函 数 、Wilcoxon 秩 和 检验 和 工程 优化 问题 上 来 测试 所 改进 算法 的 寻 优 能 力 。 实 
25s 相 比 其 他 新 提出 的 智能 算法 ，DEtanAO 算法 具有 更 强 的 寻 优 能 力 和 更 快 的 收敛 速度 。 
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Aquila optimizer integrating differential mutation and tangent flight 


Xu Yifeng, Liu Sheng, Liu Yusong, Zhang Weikang 
(School of Management, Shanghai University of Engineering Science, Shanghai 201620, China) 


Abstract: The Aquila optimizer (AO) , although capable of robust global exploration, had problems with inadequate local 
development. In the study, a differential evolution mutation and tangent flight Aquila optimizer (DEtanAO) was proposed. 
Firstly, the mutation operation in the differential evolution algorithm can make the algorithm have strong development ability 
and make up for the shortcomings of the AO algorithm. Then, the tangent flight strategy in the tangent search algorithm has a 
strong ability to explore the search space and can make the algorithm jump out of the local optimal solution, which is used to 
replace the Levy flight in the AO algorithm. The combination of these two strategies effectively balanced the exploration and 
exploitation stage of the DEtanAO algorithm. Finally, in order to verify the optimization performance of the DEtanAO 
algorithm, the optimization ability of the improved algorithm was tested in 12 standard benchmark functions, high- 
dimensional functions, Wilcoxon rank-sum test, and engineering optimization problems. The results showed that the 


E" optimization ability and convergence speed of the DEtanAO algorithm was better, compared with other newly proposed 
r intelligent algorithms. 
Key words: aquila optimizer; differential evolution mutation; tangent flight; exploration and exploitation 


®© o 引言 大 的 全 局 勘探 能 力 、 搜 索 效率 高 、 收 敛 速度 快 等 优点 ， 但 其 
p= m 局 部 开发 能 力 不 足 ， 容 易 陷 入 局 部 最 优 。 王 爽 等 四 提出 了 一 
= 随 着 现实 世界 中 离散 、 无 约束 性 质 的 优化 问题 越 来 越 复 。 种 改进 的 混合 天 鹰 优 化 器 和 哈里 斯 座 眼 优化 器 IHAOHHO)， 


杂 ， 传 统 的 基于 数学 的 优化 算法 很 难 找 到 最 优 解 ， 并 且 存 在 。 ”将 HHO 算法 的 开发 策略 集成 到 AO 算法 的 探索 策略 中 ， 并 
一 定 的 缺点 和 局 限 性 ， 例 如 收敛 至 局 部 最 优 和 未 知 的 搜索 空 。” 添加 反 向 学 习 策略 (ROBL) 以 避免 局 部 最 优 停滞 。 同 时 ， 非 
间 并 且 对 于 优化 问题 的 解决 方案 单一 等 。 为 解决 上 述 问题 ， 线性 逃逸 能 量 参数 平衡 了 算法 的 探索 和 开发 阶段 。 这 些 改进 
近年 来 ， 基 于 仿生 学 的 元 启发 式 优化 算法 不 断 被 提出 和 改进 ， ”使 得 IHAOHHO 算法 的 寻 优 性 能 得 到 提升 ， 除 此 之 外 ， 
并 以 其 操作 简便 、 参 数 调整 简单 等 特点 大 受 学 者 青睐 ， 如 : Shubham Mahajan 等 人 将 AO 算法 和 算数 优化 算法 相 结合 
加 入 自 适应 权重 模拟 鸟 群 疯 食 行为 的 混沌 粒子 群 算法 (AO-AOA) 来 解决 优化 问题 ， 并 在 23 个 基准 函数 上 验证 了 所 


ÜlChaotic-Particle Swarm Optimization, CPSO), SRE Wi GI 提 算 法 的 有 效 性 与 优越 性 ， 王 世 成 外 等 人 提出 了 一 种 改进 的 
有 反哺 行为 启发 提出 的 麻雀 搜索 算法 3(Sparrow Search Algorithm, 多 目标 天 鹰 优化 器 ， 并 用 其 设计 了 一 款 基 于 混合 固体 氧化 物 
SSA) ， 受 海洋 适 者 生存 理论 启发 提出 的 海洋 捕食 者 算法 燃料 电池 (SOFC) 系 统 的 混合 动力 系统 ， 从 而 改善 了 整体 成 
BJ(Marine Predators Algorithm,MPA)， 融 合 了 精英 反 向 和 黄金 正 本 和 用 火 效率 ; Eric Ofori-Ntow Jnrg 等 人 将 离散 小 波 变换 
弦 策 略 的 鲸鱼 优化 算法 所 (Elite Opposition Based Golden-Sine ”(DWT)、 混 沌 理论 的 相 空 间 重 构 (PSR)、 天 认 优 化 算法 (AOA) 


li: 


Whale Optimization Algorithm ，EGolden-SWOAI) 等 等 。 这 些 算 和 反 向 传播 神经 网 络 (BPNN) 混 合 起 来 提出 了 DWT-PSR- 
法 的 提出 为 解决 复杂 的 优化 问题 提供 了 更 多 的 思路 。 AOA-BPNN 模型 ， 并 将 其 应 用 于 预测 风速 ， 实 验 结果 表明 

天 订 优 化 器 (Aquila Optimizer，AO)G5 是 Laith Abualigah 该 模型 预测 精度 较 高 。 但 是 由 于 天 认 优 化 器 提出 的 时 间 较 短 ， 
等 人 于 2021 年 提出 了 一 种 新 的 基于 种 群 的 优化 方法 ， 该 方 到 内 并 没有 对 AO 算法 的 改进 进行 充分 研究 ， 其 性 能 仍 有 i 
法 受到 自然 界 中 天 鹰 在 捕获 猎物 过 程 中 行为 的 启发 ， 具 有 强 一 步 的 提升 空间 以 便 应 用 于 更 多 的 实际 问题 中 。 本 文 提 出 一 
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录用 定稿 REA, F: 


种 融合 差分 变异 和 切线 飞行 的 天 座 优 化 器 (DEtanAO)， 首 先 
在 AO 算法 中 加 入 差分 变异 策略 ， 提 高 算法 的 搜索 精度 ， 接 
着 利用 切线 搜索 算法 00(Tangent Search Algorithm ，TSA) 中 
的 切线 飞行 策略 替代 AO 算法 中 的 莱 维 飞行 策略 ， 进 一 步 优 
化 AO 的 寻 优 性 能 ， 使 其 跳出 局 部 最 优 解 。 为 证 明 改 进 算法 
的 有 效 性 ， 随 后 在 单 模 态 、 多 模 态 和 高 维 测试 函数 上 进行 多 
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| rQ+p)xsin[ 至 ] 
r( EE) o2) 
个 搜索 方法 (X,) 生成 的 + 的 下 一 
lew(D) 是 莱 维 飞行 分 布 函数 ， 
s 为 固定 为 1.5 的 


中 ， X,(1+1) 是 | Zu 
ERRE D 表示 维度 空间 ， 
X,(r) 是 取 值 范围 在 [1,N] 的 随机 解 。 其 中 ， 


种 算法 的 对 比 实验 ， 并 进行 Wileoxon 秩 和 检验 分 析 ， 进 


步 在 工程 应 用 问题 上 检验 所 提 算 法 的 前 沿 性 ， 结 果 表 明 
DEtanAO 算法 的 寻 优 精度 、 收 敛 速度 有 着 显著 提升 。 


1 “天 鹰 优 化 器 


天 鹰 优化 器 模拟 了 天 鹰 对 不 同 猎物 的 不 同 捕猎 方式 。 天 
鹰 对 于 快速 移动 猎物 的 狩猎 方式 反映 了 算法 的 全 局 探索 能 力 ， 
对 于 慢 速 移动 猎物 的 狩猎 方式 反映 了 算法 的 局 部 开发 能 
AO 算法 具有 较 强 的 全 局 探索 能 力 、 较 高 的 搜索 效率 和 较 忆 
的 收敛 速度 ， 但 其 局 部 开发 能 力 不 足 ， 容 易 陷入 局 部 最 优 。 
AO 算法 模拟 了 天 鹰 在 狩猎 期 间 的 行为 ， a 
方法 表示 : MLE AS EN AAE RRR EN (X); G8 
过 短 滑翔 攻击 的 等 高 飞行 在 搜索 空间 内 探索 eh ; 通过 慢 
速 下 降 攻击 的 低空 飞行 在 收敛 搜索 空间 内 探索 (X) ， 以 及 
通过 行走 和 抓 取 猎物 进行 俯冲 (X4) 。 
1.1 拓展 探索 OG) 
在 第 一 种 方法 (X,) 中 ， 天 座 识 别 猎物 区 域 并 通 坟 
弯 腰 的 高 拓 闭 来 选择 最 佳 狩猎 区 域 。 在 这 里 ， 天 | Eo 
翔 来 确定 搜索 空间 的 区 域 即 搜索 猎物 在 哪里 。 
ERA AAE AS ETA. 


H 
mi 


图 1 XIEIGUBXURISE ECTS EKTA 
Fig. 1 The behavior of the Aquila high soar with the vertical stoop 


这 种 行为 用 数学 式 (1) 表 示 为 
X(t+1) = Xs (»(1-z)«at. (t) - Xy (0)) rand (1) 


33,0) i 21 Dim Q) 


其 中 ， 和 4+D 是 由 第 一 个 搜索 方法 OG) 生成 的 { 的 下 一 次 迭 
UE. Xu (0) 是 第 t 次 迭代 的 最 优 解 ， 这 反映 了 猎物 的 大 
致 位 置 。(1-w7) 表示 通过 迭代 次 数控 制 探索 ，Xw 0) 表示 当 


Xu ()= 方 


常量 值 ，y 和 x* 在 搜索 中 呈现 螺旋 形式 ， 其 计算 公式 如 下 : 
y=rxcos(0) (6) 
x-rxsin(0) (7) 
r=n+U x D, (8) 
0-—-0x D,4*0, (9) 

其 中 ， 0 =(3xz)/2 ， ”为 介 于 1 和 20 之 间 的 固定 周期 指数 ， 


U 的 值 为 0.00565，D, 是 1 到 搜索 空间 长 度 的 整数 ，o% 的 值 
为 0.005。 
Levy(D) 
» 
X(t- 1) l X(t) 
: Xr(t) TIT 
t 
Xbest(b 
图 2 ”天 鹰 短 滑 翔 攻击 的 等 高 飞行 的 行为 


Fig.2 The behavior of the Aquila contour flight with short glide attack 
1.3 扩大 开发 OG) 

在 第 三 种 方法 (X,) 中 ， 当 天 座 锁 定 了 捕食 区 域 ， 天 座 
准备 好 着 陆 和 攻击 ， 随 后 垂直 下 降 并 进行 初步 攻击 来 试探 猎 
物 反映 ， 这 种 行为 称 为 低空 飞行 和 慢 速 下 降 攻击 。 图 3 展示 
了 天 磨 低空 飞行 和 缓慢 下 降 的 攻击 行为 。 


Xbest (t) 


图 3 KERE CEHROHS FEIER: GO 
Fig. 3 The behavior of the Aquila low flight with slow descent attack 


前 解决 方案 在 第 t 次 迭代 时 的 平均 值 。rand 是 [0,1] 之 间 的 随 
机 数 。t 和 了 T 分 别 表示 当前 迭代 次 数 和 最 大 友 代 次 数 。Dim 
是 问题 的 维度 ，N 是 候选 解 的 个 数 。 

1.2 缩小 探索 范围 (X;) 

在 第 二 种 方法 (X,) 中 ， 当 天 座 从 高 空 找到 猎物 区 域 时 ， 
会 在 目标 猎物 上 方 盘旋 ， 准 备 好 发 动 攻 击 。 这 种 方法 称 为 短 
滑翔 攻击 的 等 高 飞行 。 此 时 ， 天 鹰 优化 器 狭窄 地 探索 目标 猫 
物 的 选 定 区 域 ， 为 攻击 做 准备 。 图 2 展示 了 带 有 短 滑翔 攻击 
的 天 鹰 等 高 飞行 的 行为 。 这 种 行为 用 数学 式 (3) 表 示 为 

X; (t1) 2 X4, (t)xlevy(D)* X5 (t) - (y x)xrand (3) 


levy(D)zsx x sa 
viz 


(4) 


这 种 行为 用 数学 式 (10) 表 示 为 
X. (t1) (X4, (1) - Xy (t)) xa — rand 
+((UB — LB)x rand + LB)xó (10) 
其 中 ， (+D 是 由 第 三 个 搜索 方法 (X;) 生成 的 t 的 下 一 
迭代 的 解 ， 和 5 是 开采 调整 参数 ， 其 值 较 小 ， 在 (0,1) 的 范 
HA, LB, UB 分 别 代表 所 给 问题 的 上 下 限 。 
1.4 缩小 开发 范围 (X) 
在 第 四 种 方法 (X) 中 ， 当 天 鹰 接近 猎物 时 ， 它 会 根据 


猎物 的 随机 移动 攻击 猎物 ， 这 种 方法 称 为 “行走 joa i 
如 图 4 Bras. 
这 种 行为 用 数学 式 (11) 表 示 为 
X, (t1) 2 QF x X pa (t) -(G, x X (t)xrand)—G, xlevy(D) (11) 


2xrand -l 
QF (1) =1 F (12) 


G,-2xrand-l (13) 
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G-2x[ 1 二] (14) 


其 中 ， X00 是 由 第 四 个 搜索 方法 (&) 生成 的 t 的 下 一 次 
迭代 的 解 ，QEF 表示 用 于 平衡 搜索 策略 的 质量 函数 ， 其 中 
Fe(0,1) G 表示 猎物 在 逃逸 过 程 中 天 座 采 用 的 不 同方 法 ， 
GG 表示 从 2 到 0 的 递减 值 ， 代 表 了 天 鹰 在 追踪 猎物 时 从 第 一 
个 位 置 到 最 后 一 个 位 置 时 的 飞行 坡度 。 OF(D) 是 第 t 次 迭代 时 
的 质量 函数 值 ，rand 是 (0,1) 的 随机 数 ，t 和 了 分 别 表 示 当 前 
和 最 大 的 迭代 次 数 。 图 5 是 AO 算法 流程 。 


X(t) 


X(t/- 1) 


A 
w 
N * j 
LEE a 
Xpest (0 
图 4 行走 并 抓 住 猫 物 进 行 俯冲 的 行为 
Fig.4 Behavior of the Aquila walk and grab prey 


eH 
em 否 
ii | 


采 索 阶段 开发 阶段 
E 选 公 式 (10) 或 (11) 
| 


D 
d 
m 
M 
g 


法 流程 

Fig.5 The flow chart of AO algorithm 
在 AO 算法 的 两 个 阶段 ， 都 需要 添加 贪心 算法 来 保证 个 
体 的 新 位 置 优 于 其 原 位 置 。 从 AO 算法 的 实现 公式 来 看 ， 
个 位 置 更 新 公式 中 ， 式 (1) 和 (11) 会 向 0 ki, RAOS 
到 常数 ， 只 有 式 (3) 勉 强 会 收敛 至 最 优 解 ， 故 AO 算法 易 陷入 
局 部 最 优 。 


2 ”改进 的 天 应 优化 器 (DEtanAO) 


2.1 差分 变异 策略 
差分 进化 算法 UU(DE) 是 一 种 基于 种 群 进化 的 算法 。 通 
过 对 种 群 采取 三 种 进化 操作 : 变异 操作 、 交 又 操 作 以 及 选择 
操作 进行 迭代 使 算法 趋 于 全 局 最 优 解 。DE 算法 利用 差分 这 
种 变异 操作 来 产生 变异 种 群 作为 新 的 个 体 ; 
对 变异 种 群 和 原始 种 群 进行 交叉 ， 得 到 交叉 种 群 ， 对 原始 种 
群 和 交叉 种 群 利用 贪 禁 选择 选取 下 一 代 种 群 。 差 分 变异 指 的 
是 DE 通过 两 个 不 同 父 代 向 量 做 差 ， 实 现 变异 个 体 与 差分 向 
量 的 结合 ， 常 用 的 变异 策略 如 表 1 所 示 。 
A 1 中 Ps(0D HERAT, nzncznBHn.n.5 
为 随机 整数 ，Xiwj 表示 第 j RERNE, Vi 表示 变异 后 的 
个 体 ，Xi; 表示 第 7 REISDE, XL. XL. XL 表示 当前 
种 群 中 随机 选取 的 不 同 个 体 。 
IEF “best” HERAN 


TE 


有 良好 的 开发 能 力 和 快速 收 


敛 的 性 能 ， 以 天 鹰 当前 种 群 中 适应 度 最 好 的 个 体 作为 基 矢 量 ， 


RPX, X, 


天 鹰 种 群 随机 产生 。 其 中 ， 变 异 算 子 F 随 
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着 迭代 次 数 的 增加 动态 随机 变化 ， 在 算法 迭代 初期 ， 玉 值 较 
小 时 ， 缩 放 因子 对 变异 的 影响 较 小 ， 天 座 的 变异 受 群 体 中 最 
优 天 座 的 影响 ， 变 异 的 天 认 会 快速 朝向 群体 最 优 个 体 靠近 ; 
在 算法 运行 后 期 , F 值 较 大 时 ， 缩 放 因子 对 变异 影响 较 大 ， 
天 认 聚 集 在 最 优 天 鹰 个 体 的 周围 ， 使 得 种 群 在 最 优 个 体 周围 
进行 局 部 细致 探索 ， 提 高 算法 的 搜索 精度 。 
表 1 变异 策略 
Tab.1 Variation strategy 


变异 策略 数学 表达 式 
DE/rand/1 Vi =X; F(XL 7X1) 
DE/best/1 Vij = Xi; +F (XE, 一) 


DE/current-to-best/1 Vi; = Xi +F (Xis -X()«F(X;; - 


再 通过 交叉 操作 ， 


2.20 切线 飞行 策略 

切线 搜索 算法 00(The Tangent Search Algorithm，TSA) 是 
2021 年 提出 的 一 种 新 的 基于 种 群 的 优化 算法 。TSA 使 用 
个 基于 切线 函数 的 数学 模型 将 给 定 的 解 移 向 更 好 的 解 。 无 论 
是 基于 导数 还 是 基于 微分 的 最 优化 算法 都 是 类 似 于 下 式 (15) 
的 下 降 等 式 : 


X“ 2 X! estep*d (15) 
其 中 step 是 移动 的 大 小 ，d 是 移动 的 方向 。 最 优化 算法 的 区 
别 在 于 如 何 计算 步 长 ， 基 于 导数 的 方法 使 用 梯度 下 降 法 或 海 
森 和 矩阵 ， 然 而 基于 微分 的 方法 中 ， 例 如 元 启发 式 算法 ， 它 们 
使 用 随机 步骤 收敛 到 全 局 最 优 。 比 如 ， 遗 传 算法 0 使 用 高 
斯 变异 作为 步 长 ， 差分 进化 算法 0 通过 当前 种 群 个 体 之 间 
的 差 值 来 计算 步 长 ， 布 谷 鸟 搜 索 算 法 [是 通过 莱 维 飞行 函 
数 来 计算 步 长 。 
步 长 的 优化 对 于 算法 的 优化 极为 关键 ， 步 长 的 大 值 有 利 
于 探索 ， 小 值 有 利于 开发 。TSA 算法 中 提出 了 一 种 基于 切 
线 函数 的 新 的 步 长 ， 称 之 为 切线 飞行 ， 切 线 函数 有 助 于 有 郊 


地 探索 搜索 空间 。 全 局 和 局 部 游 走 相 结合 的 探索 搜索 方程 
下 式 (16) 所 表示 。 
X! = X! e step*tan(0) (16) 


从 数学 逻辑 上 讲 ， 越 接近 2890, ， 切 值 越 大 ， 得 到 的 
解 离 当 前 解 越 远 ， 越 接近 0 的 2 ， 切 值 越 小 ， 得 到 的 解 越 接 
近 当 前 解 。 探 索 搜 索 方 程 应 用 于 概率 为 VD 的 每 个 变量 ， 其 
中 了 为 问题 的 维 数 。 
切线 飞行 函数 替代 天 鹰 优化 器 中 的 莱 维 飞行 函数 使 
AO 在 规避 局 部 开发 不 足 的 同时 又 进一步 提高 AO 算法 的 全 
局 探索 能 力 ， 6. 7 是 用 Mantegna 方法 05 模 拟 1000 次 的 
莱 维 飞行 和 切线 飞行 的 随机 游 走 示 意图 。 

10 | 


50  -200  -150 -100 -50 0 50 
图 6 莱 维 飞行 随机 游 走 示意 图 
Fig.6 A diagram ofthe Levy Flight random walk 

从 图 6. 7 可 以 看 出 莱 维 飞行 产生 的 步 幅 随机 性 大 、 步 
长 区 间 窜 ， 导 致 传统 算法 可 能 出 现在 迭代 前 期 搜寻 距离 过 小 ， 
而 在 迭代 末期 搜寻 出 现 搜寻 距离 过 大 的 问题 ， 使 得 算法 优化 
和 迭代 周期 过 长 ， 精 度 不 足 ， 然 而 切线 飞行 出 现 大 跨 步 的 频率 
要 比 莱 维 飞 行 高 ， 弥 补 了 莱 维 飞行 搜寻 距离 过 大 过 小 的 问题 ， 
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切线 飞行 比 莱 维 飞行 更 有 利于 跳出 局 部 最 优 解 并 进行 大 范围 1.60GHz 2.11 GHz， 内 存 16.0GB， 主 频 2.11GHz， 仿 真 软 件 
的 搜索 ， 使 得 天 鹰 获得 更 多 的 捕食 机 会 。 为 MATLAB 2020b。 
2.3 DEtanAO 的 时 间 复 杂 度 分 析 3.2 比较 对 象 和 参数 设置 

时 间 复 杂 度 分 析 通 党 依赖 于 三 个 规则 :种群 初始 化 、 计 本 文选 取 基 本 天 鹰 优 化 器 (AO)5、 混 沌 粒子 群 算法 


算 适 应 度 函 数 和 更 新 解 。 假 设 种 群 数量 为 N，O(N) 是 种 群 初 “(CPSO)D、 麻 淮 搜 索 算法 (SSA)D、 海 洋 捕食 者 算法 (MPA)PI、 
始 化 过 程 的 计算 复杂 性 ，O(TXN)+O(TXNXDim) 是 解 的 更 融合 了 精英 反 向 和 黄金 正弦 策略 的 鲸鱼 优化 算法 (EGolden- 
新 过 程 的 计算 复杂 性 ， 其 中 ， 总 和 欠 代 次 数 称 为 T， 问 题 维度  SWOA) 员 以 及 用 切线 飞行 策略 蔡 换 莱 维 飞行 的 天 鹰 优化 器 
为 Dim。 因 此 ， AO NRN HEREA O(N X (TADI (tanAO)， 与 本 文 所 提 的 融合 差分 进化 和 切线 飞行 的 天 认 优 


1990 化 器 (DEtanAO) 进 行 对 比 。 为 保证 实验 的 公平 有 效 ， 所 有 算 
Di ——————— 4 VEIERREEXUE UE BL 30, ERAOJ 500, Eth ec 
置 如 表 2 所 示 。 
idi MACHEN. 2s X2 各 算法 的 实验 参数 
-2000 b ics € Tab.2 Experimental parameters of each algorithm 
3000 算法 参数 
AO &=0.1, 656=0.1, U -0.0056, «c-0.005 
-4000 Oo | a tanAO 2 —0.1, ó-0.1, U -0.0056, «c =0.005 
"" m^ ] | | DEtanAO a=0.1, à-0.1, U —0.0056, € 0.005 
CPSO C,21.5, C,-1.5, u-2 
6009 oo 0 100 200 300 400 SSA ST=0.8, PD-0.2, SD-0.2 
图 7 切线 飞行 随机 游 走 示意 图 MPA P=0.5, FADs=0.2 
Fig.7 Tangential flight random walk schematic EGolden-SWOA b=1 
标准 AO 的 复杂 度 可 知 ，DEtanAO 仅 添 加 了 差分 进化 3.3 测试 函数 
算法 中 的 变异 策略 O(TXN)， 并 将 原来 的 Levy 飞行 策略 变 为 验证 改进 算法 有 更 好 的 寻 优 性 能 ， 本 文选 取 12 个 基 
成 tangent 飞行 策略 ， 不 增加 计算 复杂 性 。 因 此 ，DEtanAO 准 测 试 函数 进行 函数 优化 对 比 实验 ， 其 中 选取 单 峰 测试 基准 
的 总 的 时 间 复 杂 度 为 O(NX(TX(D+1)+1))。 函数 ( ~ 所 ) 用 于 检查 算法 的 局 部 开发 能 力 ， 多 峰 测 试 基 
综 上 分 析 ，DEtanAO 算法 复杂 度 与 标准 AO 算法 是 时 间 函数 (万 ~ 记 ) 用 于 检查 算法 的 全 局 勘探 能 力 ， 固 定 维 度 多 峰 
复杂 度 属 于 同一 级 别 ， 并 未 增加 额外 的 计算 阶级 。DEtanAO RU fo~ fi;)， 具 体 函 数 信 息 见 表 3。 
算法 的 伪 代 码 如 下 所 示 。 表 3 ”基准 测试 函数 
算法 Pseudo-code of the DetanAO Tab.3 Benchmark function 
输入 : i,e,Q,6,etc 编号 函数 维 数 范围 最 优 值 
输出 : best position、best fitness Sphere dim [-100,100] 0 
Initialization the population % 种 群 初始 化 Schwefel 2.22 dim [-30,30] 0 
Define initial parameter( i,e,Q,6,etc ) % 定 义 初始 参数 f Schwefel 1.2 dim [-100,100] 0 
while(The end condition is not met) do fa Schwefel 2.21 dim [-100,100] f 0 
Calculate the fitness function values and the best Í Rosenbrock dim [-30,30] 0 
obtained solution X,(r) . % 计 算 最 佳 适应 度 值 f Quartic dim [-1.28,1.28] 0 
for (i=1,2,.…, N) do J Rastrigin dim [-5.12,5.12] 0 
Update the X„(t).x,y,G,,G,,Levy(D), ett % 更 新 各 参数 "n Ackley dim [-32.32] 0 
Expanded exploration (X,) : Update the current h Penalized dim [-50,50] 0 
solution using Equation(1) % 拓 展 探索 ， 用 式 (1) 更 新 位 置 1 fi Foxholes 2 [-65,65] ] 
Narrowed exploration (X,) : Update the current fa Branin 2 [-5,5] 0.398 
solution using Equation(3) % 缩 小 探索 范围 : 用 式 (3) 更 新 位 置 2 fia Goldstein Price 2 [-2,2] 3 
Expanded exploitation (X;) : Update the current 3.4 寻 优 精度 分 析 
solution using Equation(10) % 扩 大 开发 : 用 式 (16) 更 新 位 置 3 将 本 文 所 提出 的 DEtanAO, 5 AO. tanAO, CPSO, 
Narrowed exploitation (X,): Update the current solution SSA, MPA, EGolden-SWOA 分 别 在 12 个 基准 函数 上 独立 
using Equation(11) % 缩 小 开发 范围 : 用 式 (11) 更 新 位 置 4 运行 30 次 ， 其 中 非 国定 维度 函数 的 维度 设置 六 
end for dim-30/100/500. 3&4. 5 统计 了 7 种 算法 分 别 在 低 维 、 高 维 
for 和 国定 维 数 时 ， 得 到 的 最 优 解 、 平 均值 和 标准 差 。 
Mutation:generate a mutant population % 差 分 变异 策略 K 4 的 测试 结果 显示 ，DEtanAO ÆRA f. f. 
end for fa 上 的 寻 优 效果 达到 100%， 可 以 直接 搜索 到 对 应 函数 的 最 
end while 优 值 。 相 较 于 其 他 对 比 算 法 ，DEtanAO ERAMA 5h. ha 
function tangent flight % 切 线 飞 行 策略 万、 上 的 标准 差 为 0， 说 明 DEtanAO 具有 较 强 的 鲁 棒 性 ， 
end 差分 变异 策略 保证 了 种 群 的 多 样 性 。 在 低 维 条 件 下 ， 除 了 直 
3 ”仿真 实验 与 结果 分 析 Be * 值 的 函数 ，DEtanAO 相 比 其 他 对 比 算法 至 少 
高 出 11 个 数量 级 、 至 多 高 出 235 个 数量 级 。 在 高 维 条 件 下 ， 
3.4 仿真 实验 环境 DEtanAO 相 比 其 他 对 比 算法 在 求解 大 规模 函数 时 ， 仍 能 保 
本 次 仿真 测试 环境 为 : 操作 系统 版 本 为 Win10、64 位 操持 较 高 水 平 ， 最 大 程度 的 接近 函数 最 优 值 ， 其 中 六 和 户 仍 
作 系统 ， 处 理 器 为 Intel? Core™ i5-10210U CPU @ ”能 直接 找到 最 优 解 ， 太 保持 与 低 维 条件 下 相同 的 水 平 ， 表 


明 所 提 算 法 在 求解 大 规模 问题 时 未 
Eo (RFK S5, AEE 
以 直接 搜寻 到 最 优 
上 所 述 ，DEtanAO 


较 强 的 稳定 性 
好 ， 但 是 DEtanAO 仍 可 
看 ， 鲁 棒 性 也 较 好 。 综 J 
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陷入 “ 维 数 灾难 ”， 


Tab.4 Performance comparison on low-dimensional and 
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的 稳定 性 以 及 精 


得 算法 跳出 局 音 


RA 低 维 、 高 维 函数 上 的 性 能 比较 


EEUU. Sess $E 


high-dimensional functions 


有 度 均 较 其 他 算法 有 大 幅度 提升 。 这 说 明 ，DEtanAO 通过 引 
了 原 算法 过 度 依赖 天 鹰 根据 不 同 攻 击 方 
更 新 的 状况 ， 利 用 切线 飞行 公式 ， 使 
的 寻 优 质量 。 


、 dim=30 dim=100 dim=500 
idi 最 优 解 平均 值 标准 差 最 优 解 平均 值 标准 差 最 优 解 平均 值 标准 差 
DEtanAO 0 0 0 0 0 0 0 0 0 
tanAO 6.34E-174 1.18E-107 4.26E-108  9.01E-167  L.67E-117 9.12E-117 3.47E-163  4.70E-108 — 2.25E-107 
AO 1.32E-160  3.45E-105 1.89E-106  1.61E-157 1.15E-98 。 6.29E-98  2.33E-160  3.88E-103  2.13E-102 
f CPSO 4.98E-06 1.0769 2.7904 224E-06 0.68775 2.1899 1.22E-06 0.54379 2.0649 
SSA 0 3.15E-58 —— 1.63E-57  343E-143  445E-56 242E-55 0 528E-56  2.85E-55 
MPA 2.54E-24 — 5.38E-23  7.51E-23 8.95E-21 2.54B-19 2.25E-19  126E-17  800E-17  493E-17 
EGolden-SWOA 0 0 0 0 0 0 0 0 0 
DEtanAO 8.86E-234 — 1.50E-210 0 6.15E-231  2.19E-201 0 1.23E-228  4.87E-206 0 
tanAO 1.46E-81 — 2.61E-54 1.42E-53 1.04E-80  3.46E-32 133E-51  2.60E-78 1.96E-53  1.07E-52 
AO 3.59E-80 — 2.09E-74  449E-74 1.02E-79  324E-51  1.77E-50  L80E-78  243E-55  133E-54 
h CPSO 2.88E-07 0.28538 1.5558  4.08E-05 0.26831 1.4591 2.52E-05 0.26964 1.0202 
SSA 6.66E-177 7.29E-28 3.99E-27 2.66E-88 8.45E-26 4.30E-25 6.99E-78 7.82E-31  3.03E-30 
MPA 1.44E-14 — 245E-13  195E-13 1.74E-12  L87E-ll  L87E-ll 4.51E-11 6.19E-10 — 9.86E-10 
EGolden-SWOA  1.38E-220  8.98E-174 0 1.13E-208 —3.65E-171 0 2.59E-204 ”4.88E-161 — 2.67E-160 
DEtanAO 0 0 0 0 0 0 0 0 0 
tanAO 1.76E-159 1.17E-101 6.41E-101  3.27E-160  1.29E-101 — 7.09E-08 5.14E-101 5.01E-104 — 2.74E-103 
AO 4.96E-154 — 7.05E-103  3.67E-102  2.98E-155 2.08E-98 1.09E-100 5.20E-98 6.21E-98  3.37E-97 
h CPSO 2.41E-09 0.52855 1.7004  3.61E-05 0.90048 2.4076  603E-07 0.26886 1.018 
SSA 1.65E-299  6.72E-20  2.94E-28 0 415E-23 — 226E-22 0 5.22E-23  225E22 
MPA 2.54E-08 — 0.0001959  0.0005031  0.0018476 13.018711 36.046243 165.03471 5865.8164 4933.6315 
EGolden-SWOA 0 0 0 0 7 ATE-286 0 0 243E-273 0 
DEtanAO 2.75E-235  1.47E-216 0 145E-237  5.05E-217 0 1.25E-234  1.00E-220 0 
tanAO 2.89E-81 ^ 5.25E-54  2.87E-53 4.08E-81 ^ 4.48E-55  2.45E-54 1.69E-80 2.81E-53  L54E-52 
AO 5.25E-82 — 6.87E-55  3.76E-54 1.50E-81  9.32E-57 5.10E-56 4.57E-80 ^ 624E-82 34QE-61 
万 CPSO L17E-07 1.2247 2.8645 2.82E-06 0.51696 1.6558 544E-07 0.95575 2.5679 
SSA 0 9.95E-34 — 5.15E-33 — 5.99E-160 2.04E-32 1.11E-31 1.73E-74  1.08E31 4.32E-31 
MPA 5.24E-10 — 323E-09  147E-09  4.65E-08  229E-07  136E-07  5.83E-06  223E-05 1.32E-05 
EGolden-SWOA — 4.15E-211  221E-171 0 0 7.47E-286 0 1.17E-206 2.920E-179 0 
DEtanAO 4.34E-28  0.0001801  0.0006804 6.40E-28 ^ 2.16E-05  4.89E-05 1.75E-26  1.01E-06 ^ 2.75E-06 
tanAO 1.01E-06 — 0.0022965  0.0038784  3.81E-06 0.0016667  0.0053777  1.00E-05  0.0051573  0.0066908 
AO 1.35E-06 0.0060161  0.0108928  0.0001869  0.0232199  0.0363135  1.05E-06  0.1357954  0.2335892 
Í CPSO 420E-06 ^ 0.54473 2.05288 — 7.86E-06 — 0.27242 1.4861 1.11E-05 0.54892 1.7798 
SSA 1.33E-08 — 3.46E-05  9.84E-05 1.49E-06 0.0071157 0.0228404 6.24E-07 0.0179906 0.0352481 
MPA 24.233312 25.357126 0.5291385 95.903627 97.181244 0.6738597 496.48333 497.07268 0.1888052 
EGolden-SWOA 2.27E-11 6.69E-06 3.55E-05 2.63E-09 8.95E-07 1.89E-06 2.39E-09 8.36E-05 0.0004372 
DEtanAO 7.78E-06 0.0001963 0.0001794 8.70E-06 0.0001509 0.0001348 1.32E-05 0.0001044 6.30E-05 
tanAO 2.69E-06 7.40E-05 5.62E-05 1.10E-05  0.0001021 9.00E-05 1.12E-05 9.67E-05 7.52E-05 
AO 4.87E-06 8.34E-05 7.30E-05 4.76E-07 8.42E-05 7.89E-05 7.42E-06 7.20E-05 5.09E-05 
CPSO 9.69E-09 0.82239 2.2742  461E-06 0.27941 1.069  2.58E-05 0.56843 2.1232 
SSA 9.10E-05 —0.0017566  0.0017836  0.0001687  0.0017752  0.0012219  0.0001651  0.0033734  0.0050135 
MPA 0.0002901 0.0013413 0.0006722  0.0003923  0.0018288  0.0008126  0.0006264  0.0022338 0.0009748 
EGolden-SWOA 3.32E-06 — 8.59E-05 — 7.24E-05 5.67E-07  0.0001001 0.0001069 ^ &.48E-06 ^ 9.47E-05  7.95E-05 
DEtanAO 0 3.79E-I5 — 2.08E-14 0 0 0 0 0 0 
tanAO 0 0 0 0 0 0 0 0 0 
AO 0 0 0 0 0 0 0 0 0 
f CPSO 2.05E-06 0.54097 2.0541 L05E-08 0.2739 1.4954 &63E-08 0.54843 2.0851 
SSA 0 0 0 0 0 0 0 0 0 
MPA 0 0 0 0 0 0 0 0 0 
EGolden-SWOA 0 0 0 0 0 0 0 0 0 
DEtanAO 8.88E-16 — 8.88E-16 0 8.88E-16 — 8.88E-16 0 8.88E-16 — 8.88E-16 0 
tanAO 8.88E-l6 —— 8.88E-16 0 8.88E-l6 — 2.52E-13  1.37E-12 8.88E-16  8.88E-16 0 
f AO 8.88E-16 — 8.88E-16 0 8.88E-16 — 8.88E-16 0 8.88E-16 — 8.88E-16 0 
CPSO 2.54E-05 — 0.0017061  0.0031528 9.97E-13 0.49623 1.8854  227E-07  0.0031208 0.011286 
SSA 8.88E-16 — 8.88E-16 0 8.88E-16 — 8.88E-16 0 8.88E-16 — 8.88E-16 0 
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BRA 低 维 、 高 维 函 数 上 的 性 能 比较 
dim=30 dim=100 dim-500 
函数 算法 ui TIUS z = A " zx m 
最 优 解 平均 值 标准 差 最 优 解 平均 值 PREZ 最 优 解 平均 值 标准 闫 
MPA 3.53E-13 1.42E-12 L16E-I2  7.10E-12 4.60E-11 2.34E-11 2.22E-10 4.30E-10  1.29E-10 
' EGolden-SWOA — 8.88E-16 8.88E-16 0 8.88E-16 8.88E-16 0 8.88E-16 8.88E-16 0 
DEtanAO 2.81E-32 6.72E-08 | 2.26E-07 1.69E-32 3.17E-09 1.20E-08 1.39E-32 1.33E-10 — 5.87E-10 
tanAO 2.66E-12 7.58E-07 . 1.51E-06 | 2.53E-09 2.95E-07 4.86E-07 3.59E-10 1.89E-07 . 4.01E-07 
AO 6.56E-09 4.38E-06 | 9.48E-06 . 3.82E-09 3.07E-06 7.42E-06 1.09E-08 6.15E-07  7.44E-07 
h CPSO 1.11E-05 0.68033 2.1691 8.34E-06 0.6621 2.1033 0.0001318 0.41022 1.6491 
SSA 8.60E-18 3.6SE-12 1.81E-11 1.37E-10 7.40E-08 3.43E-07 2.49E-12 5.62E-08 1.10E-07 
MPA 1.22E-09 5.69E-09 . 3.34E-09  0.0236108 0.0441271 0.0122656 0.3516941 0.4089202 — 0.0289135 
EGolden-SWOA — 6.15E-14 4.88E-08 | 2.43E-07 . 1.24E-14 8.57E-10 3.28E-09 8.56E-14 3.11E-11 — 6.39E-11 
表 5 固定 维 数 函数 上 的 性 能 比较 
Tab.5 Performance comparison on fixed dimensional functions 
算法 fo fa f 
最 优 解 平均 值 标准 差 最 优 解 平均 值 PREZ 最 优 解 平均 值 标准 闫 
DEtanAO 0.9980038 2.6384379 2.895985  0.3978874 0.3978881 3.96E-0 3 3.0038 0.01813 
tanAO 0.9980038 . 3.7318068 4.590161  0.3979267 4 0.4022246 . 0.0131830 3.0003929 3.3083947 . 0.414067 
AO 0.9980038 . 2.7607481 3.459285  0.3978936 . 0.3980834 0.0001897 . 3.00003146 . 3.03204623 0.033251 
CPSO 3.29E-05 0.13484 0.73388 2.25E-08 0.25729 1.4014 3.29E-06 0.7202 1.8772 
SSA 0.99800384 . 3.0121167 4.260976  0.3978874 . 0.39788736 0 3 4 4.929503 
MPA 0.99800384 — 0.9980038 . 1.51E-16 0.3978874 . 0.39788736 1.23E-14 3 3 1.64E-15 
EGolden-SWOA 0.99800384 — 0.9980038 . 3.50E-12 0.3978875 0.3979736 | 0.0001675 . 3.00000002 . 3.00000887 1.11E-05 


3.5 “收敛 曲线 分 析 
算法 的 优 劣 可 以 通过 其 收敛 
敛 曲 线 展 示 了 算法 的 收敛 速度 和 陷入 局 部 最 


MPA, EGolden-SWOA 共 7 种 算法 在 30 维 情况 下 ， 对 上 述 中 可 知 ，tanAO 
12 个 基准 函数 的 收敛 曲线 对 比 图 。 观 察 上 述 图 中 各 算法 的 升 ， 观 察 收敛 曲 


线 直 观 地 表现 出 来 ， 收 


图 8(a)~Q) 列 出 了 DEtanAO、AO、tanAO、CPSO、SSA、 


这 不 仅 说 明 DEtanAO 对 
而 且 不 易 陷入 局 部 最 优 解 ， 
发 能 力 。 其 中 为 验证 改进 策略 的 
策略 替换 莱 维 飞行 策略 的 算法 记 为 tanAO， 从 收敛 曲线 
比 基 本 的 AO 算法 在 收敛 速度 上 有 明显 
线 图 可 知 tanAO 在 仅 150 次 迭代 就 达到 


优 值 的 次 数 。 


收敛 曲线 可 以 看 出 DEtanAO 的 收敛 速度 在 整 


个 迭代 过 程 中 


都 快 于 其 他 6 种 算法 ， 收 敛 精度 也 是 这 7 种 算 
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Fig.8 Convergence curve 


比 其 他 算法 有 更 好 的 全 局 探索 能 
平衡 了 全 局 探索 能 力 和 局 部 
了 效 性 ， 将 仅 用 切线 飞 和 
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了 最 优 值 ， 而 基本 AO 算法 在 250 次 迭代 才 达 到 最 优 值 ， 
tanAO 的 寻 优 精度 和 速度 比 AO 有 较 大 幅度 提升 。 
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综 上 所 述 ， 本 文 提出 的 DEtanAO 算法 不 管 在 低 维 、 高 DEtanAO 算法 的 显著 性 强 于 其 他 算法 ;， 当 p>0.05 时 ，h 显 
维 还 是 固定 维度 的 函数 下 都 具有 较 好 的 寻 优 性 能 ， 其 求解 示 为 “0”， 表 示 DEtanAO 算法 的 显著 性 弱 于 其 他 算法 ; 当 
精度 、 收 敛 速度 和 和 鲁 棒 性 均 优 于 与 之 对 比 的 其 他 六 种 算法 ， P 显示 为 “N/A” 时 ， 表 示 无 法 进行 显著 性 检验 ，DEtanAO 
有 着 明显 的 求解 优势 。 算法 显著 性 可 能 与 其 他 算法 相同 。 
3.6 Wilcoxon 秩 和 检验 将 DEtanAO 与 AO, tanAO, CPSO, SSA, MPA, 
Wilcoxon019 秩 和 检验 是 一 种 非 参 数 统计 检验 方法 ， 用 EGolden-SWOA 独立 运行 30 次 的 平均 值 进行 Wilcoxon 秩 
于 判断 比较 DEtanAO 算法 与 其 他 算法 是 否 有 显著 性 区 别 ， 和 检验 ， 从 表 6 中 的 统计 结果 可 以 看 出 ， 大 部 分 的 p 值 均 
其 结果 分 析 如 表 6 所 示 ， 其 中 “p” 表 示 检 验 结果 ,“h” 表 小 于 0.05， 表 明 DEtanAO 与 其 他 算法 存在 显著 差异 ， 并 
示 显 著 性 判断 结果 。 当 p<0.05 时 ，h 显示 为 “1”， 表 示 DEtanAO 的 性 能 在 多 个 函数 上 优 于 对 比 算法 。 
表 6 基准 函数 Wilcoxon 秩 和 检验 结果 
Tab.6 Wilcoxon rank sum test results of reference functions 
-— AO tanAO EGoldenSWOA MPA SSA CPSO 
i p h p h p h p h p h p h 
Íi 1.21E-12 1 1.21E-12 1 N/A 0 1.21E-12 1 1.21E-12 1 1.21178E-12 Į 
Í 3.02E-11 1 3.02E-11 1 1.61323E-10 1 3.02E-11 1 3.02E-11 1 3.01986E-11 1 
f 1.21E-12 1 1.21E-12 1 N/A 0 ].21E-12 1 1.21E-12 1 1.21178E-12 1 
h 3.02E-11 1 3.02E-11 1 3.01986E-11 1 3.02E-11 1 3.02E-11 1 3.01986E-11 1 
f 1.41E-09 1 1.29E-09 1 0.077271976 0 3.02E-11 1 2.13E-05 1 2.43863E-09 1 
h 0.099258 0 0.145319 0 0.016284809 1 3.69E-11 1 8.48E-09 1 5.53286E-08 1 
f 0.333711 0 0.333711 0 0.333710696 0 0.333711 0 0.333711 0 1.72025E-12 1 
h N/A 0 N/A 0 N/A 0 1.21E-12 1 N/A 0 1.21178E-12 i 
h 3.02E-11 1 3.02E-11 1 2.37591E-07 1 8.15E-11 1 9.51E-06 1 3.01797E-11 1 
fo 0.501102 0 0.841786 0 7.18381 E-05 1 5.78E-08 1 0.042542 1 1.72466E-07 1 
fü 3.16E-12 1 3.16E-12 1 3.15782E-12 1 0.024657 1 0.081493 0 3.22767E-07 1 
fo 1.36E-08 1 1.02E-10 1 1.01948E-07 1 0.186128 0 0.509236 0 1.01948E-07 1 
1/0 8/4 8/4 7/5 10/2 8/4 12/0 
3.7. 工程 优化 问题 应 用 得 算法 跳出 局 部 最 优 ， 提 高 了 算法 的 寻 优 质量 。 通 过 12 个 
为 验证 DEtanAO 算法 的 前 沿 性 以 及 在 实际 工程 中 的 优 标准 基准 函数 、 高 维 函数 、Wilcoxon 秩 和 检验 以 及 在 工程 
越 性 ， 选 取 了 压力 容器 设计 优化 问题 ， 并 通过 与 上 述 对 比 算 优化 问题 上 的 应 用 的 结果 显示 ，DEtanAO 具有 较 强 的 跳出 
法 的 改进 算法 如 融合 正 余 弦 和 柯 西 变异 的 床 省 搜索 算法 局 部 最 优 的 能 力 、 更 快 的 收 和 敛 速 度 以 及 更 高 的 收敛 精度 。 下 
(SCSSA)03、 混 合 策略 改进 的 麻 汰 搜索 算法 (MSSSA)05 基 于 一步 考 虑 将 DEtanAO 算法 应 用 到 其 他 实际 问题 的 建 模 中 ， 
分 段 式 随机 惯性 权重 和 最 优 反 馈 机 制 的 鲸鱼 优化 算法 比如 机 器 学 习 各 种 预测 分 类 算法 中 超 参 数 的 优化 ， 并 进一步 
(GEWOA)09、 高 斯 量子 粒子 群 算 法 (G-QPSO)P0 进 行 验证 比 将 优化 后 的 机 器 学 习 算法 用 来 预测 量化 投资 中 股价 变动 的 趋 
较 。 压 力 容器 设计 问题 是 一 个 经 典 工 程 优化 设计 问题 ， 势 ， 以 扩展 本 文 所 提 算 法 在 金融 等 领域 的 应 用 范围 。 
目的 是 通过 优化 圆柱 容器 身长 L、 圆 柱 容器 内 径 R BEK 表 7 压力 容器 设计 问题 求解 结果 
壁 厚 六 和 头 部 壁 厚 1 变量， 使 压力 容器 制作 成 本 达到 最 小 。 Tab. 7 Pressure vessel design problem solving results 
问题 的 目标 函数 、 约 束 条 件 和 自 变量 取 值 范围 如 下 所 示 。 算法 最 小 值 平均 值 标准 差 
和 =[D x, |HIT, T. R,L DEtanAO 5862.69 6081.35 2.1SE+02 
Minf (x) 2 0.6224x x, x, +1.7781x,x3 - 3.1661x,x? -19.84x4x2 FWOA 5886.97 - - 
g (x) 7 =x +0.0193x <0 SCSSA 8050.92 : : 
82 (x) == +0.00954x, <0 MSSSA 5897.15 6615.81 5.53E402 
S e (r)e md E 41296000 <0 (7) G-QPSO 6059.72 6342.45 4.48E+02 
g4,(x)2—x,-240x0 " 
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